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Abstract 


The development of digital controls for turbo- 
jet and turbofan engines can be facilitated by the 
use of real-time computer simulations of the en- 
gines. The engine simulation provides a "test-bed" 
for evaluating new control laws and for checking 
and "debugging" control software and hardware prior 
to engine testing, This paper describes the devel- 
opment and use of real-time, hybrid computer simu- 
lations of the Pratt & Whitney TF30-P-3 and FIOD- 
FW-100 augmented turbofans In support of a number 
of controls research programs at the Lewis Research 
Center. The role of engine simulations in solving 
the propulsion systems integration problem is also 
discussed. 


Introduction 


The development of controls for aircraft pro- 
pulsion aystcDS depends, to a great extent, on one 
being able to accurately predict the performance of 
the aircraft, the propulsion system, and their as- 
sociated controls. Computer simulations provide a 
means of analyzing the behavior and Interactions of 
these increasingly complex systems prior to full- 
scale testing of the hardware. Simulations also 
serve os aids In solving problems that arise after 
the developmdnt phase Is completed. 

Of particular interest today are simulations 
of gas turbine engines and tliclr application to the 
development of advanced propulsion system controls. 
This paper describes the development and use of 
real-time (R-T) , hybrid computer siniuloclons of the 
Pratt S Whitney TF30-P-31 and FlOO-PW-lOO^ aug- 
mented turbofans In support of a number of controls 
research programs2-5 oc the Lewis Research Center. 
In addition to the discussion of these real-time 
simulations, the role of other types of engine sim- 
ulations will also be discussed. 

Engine Simulations 


The complexity of an engine simulation is, of 
course, determined by the particular application. 

A simple, choked-oriflce rcprcscncation of a turbo- 
jet engine may be adequate when studying the dy- 
namic performance of a supersonic inlet at a given 
operatlnj condition.^ There arc, however, cases 
where an engine must be considered as more than 
Just a simple Inlet termination. For example. In 
determining the effect of engine response on an 
externally blown flap aircraft,^ u linear, trans- 
fer function representation of a two-spool turbo- 
fan engine was used. As shoi'in In Fig. 1, the 
transfer functions related components of engine 
thrust to a pilot-commanded thrust change. In this 
case, the engine simulation was combined with a 
simulation of a STOL transport In a piloted, 
ground-based simulator at the Ames Research Center. 
This transfer function approach to engine simula- 
tion can be extended to Include additional Inputs 
and outputs. Unfortunately, the adequacy of the 
transfer function models Is usually limited to op- 


eration around selected operating points due to the 
nonllnearitlcs of the propulsion system. Therefore, 
another approach must be taken when analyzing the 
gross transient behavior of engines and their eon- 
crola over a range of flight conditions end power 
settings. 

There are a number of wide-range, nonlinear 
engine simulations currently being used to support 
controls development programs. These Include both 
steady-state and transient simulations. Steady- 
state simulations are useful in establishing set- 
point control schedules end for defining engine 
limits that the control must not exceed. The most 
detailed steady-state simulations have been devel- 
oped by the various engine manufacturers for spe- 
cific engines and arc often updated to reflect 
actual engine test data. Generalized, steady-state 
programs^ are also available for simulating a 
wide range of turbojet and turbofan engine configu- 
rations. In general, steady-state simulations ara 
implemented using digital computers because of the 
digital computer's precision, repeatability, and 
f IcxlbJ llty . 

Transient simulations are also needed as design 
tools. Again, the most detailed simulations have 
been developed for specific engines and Implemented 
using digital computers. They do, however, consume 
much expensive computer time since they require 
Iterative solutions and numerical integration. 
Transient, digital simulations hsvo also been gen- 
erallzcdlO at the cost of Increased computing times. 

A similar generalized, transient pimulatlon has 
been developed for the hybrid (analog-digital) com- * 
puter.^'^ Because of its exact integration with re- 
epcct to time and compromises in modeling detail, 
this program runs faster than the all-digital simu- 
lations. However, it is not fast enough to run in 
real-time. 

Real-Time Engine Simulation 

The subject of real-time (R-T) engine simula- 
tion Is gaining increasing attention because of Its c 
application to the evaluation of new digital con- 
trol hardware and software. This concept Is 
illustrated In Fig. 2. Engine running, whether 
stocic or in a flying test-bed, is expensive and 
time-consuming. Hence, a simulation which can re- 
duce the amount of required engine testing is 
highly desirable. To serve os a realistic tost ve- 
hicle for n control, the simulation should 
(1) statically and dynamically represent the per- 
formance of all parts of the engine-control loop 
not available as hardware, (2) provide suitable 
signals at the control Interface (this requires 
digltal-to-analog, D/A ard analog-to-dlgltal , A/D 
conversions If engine simulation Is Implemented 
with a digital computer) and (3) run in real-time 
so as CO realistically Interact with the dynamics 
of the control. 

There are a number of factors to consider when 
selecting a computer for performing R-T simulations. 
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from Angular momencum equaclons. The cffucts of 
.fluid Dopicntum ov. the engine dynamics ace included 
in £)iu duct and augmentor models. 

A number of sasumptlons and model slmplifica- 
Cions ore necessary to aatisfy the R-T ccquircmcnt 
with a icusonable amount of computing equipment. 

In the TF30 end FlOO R-T slmuiationB, these incluae 
the following; (1) fan and compressor temperature 
ratios are assumed to be piecewise linear functions 
of the corresponding pressure ratios; tUlo elimi- 
nates the need for efficiency mop generation and 
exponentiation in the digital computer, (2) the 
static pressure balance between core and duct 
streams , assumed in the baseline simulations , is 
represented by a constant ratio of total pressures 
in the R-T simulations, (3) gas prop"crclea such as 
specific heats arc assumed to be constant except 
for the fan and compressor-discharge specific heats 
which arc assumed to be linear functions of temper- 
ature; the linear functions provide a good match of 
design rotor speeds along a given operating line. 

Steady-state and transient evaluations of the K-T 
simulations have shown that these assumptions do 
not significantly reduce die accuracy of tlic simu- 
lations . 

Real-Time Simula tlop Using the Hybrid Computer The equations describing the TF30 and FlOO 

"" mathematicol models have been programmed on the _ 

The hybrid computer, because it includes both Lewis Resesrcli Center's EAl Model 650 hybrid com- ^ 

a digital and analog computer, can satisfy the re- puter and EAI Model 680 analog computer. An addl- 

quireoents of R-T engine simulation. Fig. 3 shows tionai EAI Model 231R analog computer is required 
one approach to using the hybrid computer in this for the TF30 simulation because of its more compll- 

appllcar.lon. Tlie digital portion of the hybrid cated configuration. As previously stated, the 

Computer can be used to perform all of the required digital portion of the hybrid computer is used prl- 

function generation. This is the approach that is marlly for the function generstion. In each of the 

used in the TF30 and FlOO simulations discussed in R-T simulations, eight bivariate table lookups are 

the next section. BecausK a digital computer is tcqulmd atid uonsume approximately A. 7 mllllnccondo. 

used, the upper limit on update time must be ad— The function generation routine Chat is used is 

hered to. Additional calculations may be performed discussed in the next Section. In the FlOO simuln- 

digltally if these can be completed within the al~ Cion, an engine thrust calculation has been added 

lowablc update interval. The analog portion of the to the digital program resulting In a total update 

hybrid computer is available to perform the remain- time of 6.5 milliseconds. Control Inputs, such as 
ing calculations wilieh include the integration as- fuel flows and nozzle areas, are provided by digi- 

aoclated with the engine dynamics. The analog por- tal control systems programmed on the Lewis Ra- 

tion of the hybrid computer also provides the nec- search Center's SEL810B computer, 
essary control Interface since D/A and A/D con- 
versions ate performed by the hybrid computer in- HAF2 Function Routine 

terface. In addition to performing the calcula- 
tions associated with function generation, the Because of the nature of the fan and compres- 

hybrld's dlgit.il computer can also be used to auto- sot performance maps (Fig. 5) , rectilinear inter- 

mutically set-up and check-out the analog portion polatlon cannot be used to generate these functions, 

of the simulation. Therefore, a rodlal-lntetpoiation, bivariate- 

function generation routine, MAP2 has been devol- 
TF30 and FlOO R-T Simulations oped at Lewis for generating this type of function. 

The MAP2 routine is based on a routine developed 

In support of controls research programs In- eariler^^ for the same application. Tabular data 

volvlng the Pratt 6 Wltitncy TF30-P-3 and FIOO-PH- for up to twelve functions of six pairs of inde- 

100 augmented turbofans, R-T liybrld computer slmu- pendent variables are stored in a conmion block and 

latlons of those engines have been developed at shared by HAP2 and the calling program. Both 

Lewis. Fig. A illustrates, in block dlogram FORTRAN and assembly language versions of MAP2 

form, the approach taken in modeling those engines, have been developed. For the R-T simulations, the 

The engine models are patterned, as much os pos- asacmbly language version is used. An additional 

slble, after the engine manufacturer's digital entry point fWP2L has been added to the assembly 

'(baseline) simulations. The bivariate function language version to allow multiple functions of a 

data, defining the performance of the rotating single pair of independent variables to he gener- 

components, have been extracted from the digital ated with only one table search. Because ot its 

simulations. Intercoraponent volumes are assumed favorable qualities, MAP2-HAP2L is also used to 

at engine locations where either (1) gas dynamics parform the turbine function generation although 
are considered to be Important or (2) gas dynamics rectilinear Interpolation routines could be used 
are required to eliminate the need for iterative Instead, 

solutions. Time-dependent forms of the continuity, 
energy, and state equations are solved in each In- 
tcrcomponcnt volume. The rotor speeds are computed 
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If a digital computer is usud, there is an uppur 
limit on the digital update time (l.c., time re- 
quired to perform the digital calculations) that 
must not be exceeded if the closed-loop system's 
dynamic accuracy and stability are to be main- 
tained over the frequency rongo of interest. For 
most turbo fan engine simulations, this requires 
digital update times under 10 miillsecondB.li^*^2 
Achieving these update times without sacrificing 
too much in aimulation accuracy is a formidable 
toak. In addition, a general-purpose digital com- 
puter may not have the necessary interface (A/D 
and D/A) capability. 

Although they provide the necessary computing 
speed and interface, analog computers arc uaually 
not considered for turbofan engine simulations be- 
cause of the large number of multlvatUate func- 
tions Chat have to be generated. Tlicse functions 
are needed to describe the overall performance of 
the engine's rotating components (i.c., fans, com- 
pressors, turbines). The amount of required ana- 
log equipment is usually prohibitively large and 
would require considerable set-up and check-out 
time. The digital computer is, of course, well 
suited for the task of function generation. 
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R-l* Simu lat ion Verification 
I Steady-State ; 

The atcady-Btate accuracy of the R-T slmula- 
tlcna can be evaluated by running the slamlatlona 
with fixed control Inputs corresponding to the 
bosellna simulation values at selected operating 
conditions . this open-loop approach allows simu- 
lation end control Implcmeiicatlon errors to be 
Isolated. The resulting values of simulated engine 
variables such as rotor speeds can be compared with 
baseline simulation and/or engine test data. 

Fig, 6 shows results from a comparison of FlOO R-T 
simulation and experimental data at sea-level, 
static conditions. These results, together with 
results from more comprehensive comparisons,^*^ 
Indicate that the R-T hybrid simulations do ade- 
quately represent the steady-state performance of 
the TF3Q-P-3 and FlOO-PW-lOO turbofans. 

Transient 

The R-T simulations should olso accurately 
predict Che transient performance of the engines. 
One factor that can cause dynamic errors and, ulti- 
mately, Instabilities In a hybrid simulation is the 
digital update time wliich Includes the time re- 
quired to sample the analog Input variables , per- 
form the necessary calculations, and transfer the 
digital outputs to the analog computer. Tiic digi- 
tal update time appears as a time delay to the ana- 
, log computer and will generate dynamic errors If It 
is too long. This problem has been analyzed^^*^^ 
for the simple one-loop case. Unfortunately, the 
; turbofan engine presents a more complicated prob- 
lem for analysis. In the earliest R-T simula- 
tions,^'^ all analog inputs were sampled at the be- 
ginning of the digital cycle and all outputs were 
transferred to the analog after all calculatlona 
were completed. From a time delay, or phase shift, 
point of view this Is Che worst approach. 

Fig. 7 shows how the time delay can affect Che 
simulated response of Che FlOO compressor rotor 
speed to sinusoidal oscillations In main burner 
fuel flow. The effect of the digital process can 
be seen by comparing the real-time rasulcs with 
data obtained by running the simulation slower 
chan real-time. Slow-tlmc operation effectively 
speeds up the digital computer, thus minimizing Its 
effect on the frequency response. TItc comparison 
indicates that the digital contributes over 80 de- 
grees of phase shift ft 8 hertz when running in 
real-time. 

One approach to minimizing the unwanted phase 
shift is to sample analog inputs as needed and out- 
put digital data Co the analog as available. This 
approach Is Illustrated in Fig, 8. While the total 
.digital updote time Is slightly increased due to 
less efficient sampling and outputting of d.'ita, the 
reduction In calculation time for each loop can re- 
sult in a significant reductioa In the phase shift. 
For the response shuwu in Fig. 7, approximately 
one-half of the unwanted phase shift was eliminated 
by using this approach. Based on these results, 
the earlier R-T simulations have been modified 
accordingly. 

As in the case of the steady-state evaluation, 
the transient evaluation of the R-T simulations can 
best be accomplished in an open-loop fashion. The 
control inputs to the engine simulation are sched- 


luled as functions of time to match cloaed-loop, 
ibasellne simulation or experimental data. Fig. 9 
jahows a comparison of R-T simulation and expcrl- 
imental responses of FlOO engine variables to an 
Idle-to-lntermedlate power lever "slam." The ex- 
parlmental data were obtained in a Lewis altitude 
test facility which almulated the 30 000 foot , 

0.7 Hn condition. Experience with both the TF30 
and FlOO R-T almulntiona has Indicated that the 
main burner specific heat ratio has to be decreased 
by a factor of 20 to 25 to match baseline simula- 
tion reaponses. This appears to be a simple way of 
accounting for the transient effects of factors 
such as heat transfer to the engine mctal^' without 
further complicating the model. Results such as 
chose shown In Fig. 9 Indicate that this reduction 
docs result In good agreement between the R-T simu- 
lation and actual engine responses. The validity 
of this approach should be further tested as more 
experimental data becomes available, 

TP30 Digital Control Evaluation 

The R-T simulation of the TF30-P-3 engine has 
been used to evaluate a digital computer implemen- 
CaClon of the standard blll-of-materials (BOM) con- 
trol modes. ^ This evaluation Is intended to demon- 
acrace an efficient utilization of the digital con- 
trol computer' s^^ core capacity and computing time 
so as to provide sufficient capacity for extended 
control functions such as Integrated inlet and 
engine control,^® and fall-operational control.^ 
Control accuracy and dynamics , comparable to the 
hydromcchanlcol BOM system, must also be demon- 
strated. 

In partlcul.tr, rh«. R-r simulation provides a 
means of predicting (1) the stability of internal 
dynamic loops in the control, (2) Che effects of 
control calculation and update times on the engine 
performance and (3) the aLlllcy of the digital con- 
trol to perform statically and dynamically over a 
wide range of flight conditions and power settings. 

Tile R-T evaluation of the TF30 digital control 
has shown chat all of the stated objectives can be 
satisfied by the digital version of the BOM control 
modes, ^ fig. 19 shows a typical res . it obtained 
during chat cvaluaclon. The effect of increasing 
Che control update time on the response of the sim- 
ulated engine variables to an Idle-to-lr.tennedlace 
power lever "slam" is shown. For this case, the 
control update time and calculation time were equal. 
Performance degradation is observed for update 
times greater than 20 mllllseci.nds. 

Studies were also performed with a fixed con- 
trol calculation time of 3.2 milliseconds and In- 
creased update times. The idle time (difference 
between the calculation and update times) can be 
used for the previously mentioned control functions. 
Results from Chose studies Indicate that control 
update times as great as 30 milliseconds can be 
tolerated . 

Sensor Fail-Operntlonal Control Development 

To dace, Che most extensive use of the R-T 
hybrid simulations at Lewis has been In support of 
sensor fall-operational control studies.^ For 
digital control Co be viable, steps must be taken 
to ensure that the reliability of full-authority, 
flight -eiialillc'l digital control systemB matches 
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chat of the well-proven hydroncchenlcal uysteDB, 

This applies to all aspects of the control such as 
sensors, actuators, and the computer Itself. Re- 
dundant sensors and actuators may not be practical 
In certain applications due to weight and Installa- 
tion restrictions. An alternate approach for han- 
dling sensor failures la fall-operational con- 
trol^ *5 uhcrc estimates of folled-sensor outputo 
ore used by the control when failure of one or more 
sensors is detected. These estimates are based on 
Information scored In the digital control computer. 

In lieu of a full-scale engine test program, 
the R-T simulation of the TF30-P-3 engine is being 
used Co evaluate a prograsslon of fail-operational 
control schemes. Fig. 11 illustrates how a sensor 
fail-operational control works. The basic scheme 
involves estimating a failed-sensor output from 
algebraic functions of the working sensor outputs. 
The functional data can be obtained ot various op- 
erating conditions, both steady-state and tran- 
sient, and averaged appropriately. This approach 
has been integrated with the digital version of Che 
TF30-P-3 DOM control and has been evaluated on the 
R-T simulation.^ The fail-operational control adds 
a maximum of 13 milliseconda to the 3.2 millisec- 
onds required by the BOM control. Tlie fall- 
operational control (Including scored data) added 
about 6.6 K to the 5 K of core storage required by 
the BOM control. The evaluation of the fall- 
opcratlonel control using the R-T simulation has 
shown that! (1) for nonaugmcnccd operation, three 
out of the four engine sensors can be failed with- 
out violating engine limits or significantly de- 
grading engine performance, (2) for augmented oper- 
ation, only one sensor can be foiled at a time and 
a power lever rate limit has to be imposed; this 
' requirement can be attributed to the existence of a 
sensitive, integrator loop in the BOM augmenCor 
’control, and (3) satisfactory transient results are 
obtained even when the scored funcclonoA data has 
been obtained at steady-state conditions. 

Fig, 12 shows the effect of e simulated com- 
pressor rotor speed sensor failure on the response 
of engine thrust to an idle-co-maxlmum power lever 
"slam," In this case, the stored engine data was 
"learned" at a different flight condition than chat 
at which the failure occurred. The thrust response 
is somewhat slower with the failed sensor, particu- 
larly during the augmencor transient. The final 
thrust value is somewhat higher with the speed 
sensor failure than without. The augmencor fuel 
flows appear to be quite sensitive to the estimated 
speed schedule. 

Multivariable Controls Development 

Over the pest several years, aircraft opera- 
tional requlremcius have dictated the development 
of gas turbine engines with increased performance 
over a wider operating envelope. These development 
efforts have resulted in today's complex augmented 
curbofans and will result in even more complex 
variable-cycle engines in the future. Control of 
these engines will be more difficult and existing 
control design techniques, uhleh utilize loop-by- 
loop logic development, will not be able to solve 
the multivariable controls problem in an orderly 
and systematic manner. 

One opproach to the multivariable control ayn- 
Chesls problem is the application of modern control 


theory. The Air Force and NASA are currently fund- 
ing a multivariable controls development program 
for the FlOO engine. Systems Control, Inc., with 
the support of Pratt A Whitney, is applying linear 
quadratic regulotor (LQR) theory^” Co the FlOO con- 
trols problem. The resulting control will be eval- 
uated on the R-T engine simulation prior to full- 
scale engine testing. 

The LQR design procedure will utilize linear 
models of the FlOO engine. Each of the linear 
models will define the dynamic behavior of the 
engine at a specified flight condition and power 
setting. These models are being generated by Pratt 
& Whitney from basuline simulation data. Since the 
multivariable control Is to be evaluated on the R-T 
simulation, agreement beeween baseline digital and 
R-T hybrid simulations is necessary if delays In 
the program are to be avoided. One way of compar- 
ing the simulations is to generate linear models 
with each simulation and then compare transient 
responses obtained with the linear models. Pro- 
grams have been developed for automatically gener- 
ating the linear msdel coefficients from the 
sixteenth-order hybrid simulation. Fig. 13 illus- 
trates how well the resulting lineor models repre- 
sent the nonlinear simularion. The responses of 
selected engine variables to a step Increase in 
main burner fuel flow are shown. These results In- 
dicate that the linear models, generated by the 
hybrid simulation, can be used in an off-line com- 
parlsor with linear models generated by Che base- 
line digital simulation. Any differences between 
Che digital and hybrid simulations can Chen be re- 
solved prior to the multivariable control evalua- 
tion. 

Conclusiona 

Computer simulations of gas turbine engines 
play an important role in the development of ad- 
vanced propulsion system controls. In particular, 
a real-time simdlation of e turbojet or Curbofon 
engine can serve as a cesc-bed for evaluating 
actual control hardware and software prior to full- 
scale engine testing. Accurate, real-time, non- 
linear simulations are now possible due to the 
availability of modem hybrid computers. Real-time, 
hybrid computer simulations of the Pratt 6 Whitney 
TF30-P-3 and FlOO-PW-lOO augmented curbofans are 
currently being used to support a number of controls 
research programs at the Lewis Research Center. 
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Figure Z - Digital control evaluation u$ing a real-time engine simulation. 
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Figure 4. > Block diagram of FIDO mathematical model. 
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(b) MINIMUM PHASE-SHin APPROACH. 

Figure 8. - Approach to minimizing digital effects in hybrid simuiations. 
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Figure 13. - Verification of FlOO linear model at sea-level, 
static-intermediate power. 2 Percent step increase 
in fuel flow. 
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